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ABSTRACT 

While astrochemical models are successful in reproducing many of the observed in¬ 
terstellar species, they have been struggling to explain the observed abundances of 
complex organic molecules. Current models tend to privilege grain surface over gas 
phase chemistry in their formation. One key assumption of those models is that rad¬ 
icals trapped in the grain mantles gain mobility and react on lukewarm (^ 30 K) 
dust grains. Thus, the recent detections of methyl formate (MF) and dimethyl ether 
(DME) in cold objects represent a challenge and may clarify the respective role of 
grain surface and gas phase chemistry. We propose here a new model to form DME 
and ME with gas phase reactions in cold environments, where DME is the precursor 
of ME via an efficient reaction overlooked by previous models. Furthermore, methoxy, 
a precursor of DME, is also synthetized in the gas phase from methanol, which is 
desorbed by a non-thermal process from the ices. Our new modelreproduces fairy well 
the observations towards L1544. It also explains, in a natural way, the observed corre¬ 
lation between DME and ME. We conclude that gas phase reactions are major actors 
in the formation of ME, DME and methoxy in cold gas. This challenges the exclusive 
role of grain-surface chemistry and favours a combined grain-gas chemistry. 
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1 INTRODUCTION 

Relatively complex (i.e. containing more than six atoms) 
organic molecules, referred in the literature as COMs, were 
first discovered in the warm ( > 100 K) and dense ( > 10® 
cm“®) regions surrounding high mass protostars, called hot 
cores (Blake et al. 1987). Several models were then devel¬ 
oped to explain the observed abundances of those COMs 
(Millar et al. 1991; Charnley et al. 1992; Caselli et al. 1993). 
The basic idea of those models was that simple hydrogenated 
molecules like H 2 CO, CH 3 OH and NH 3 are formed during 
the cold pre-collapse phase on the cold ( < 20 K) grain sur¬ 
faces (e.g. Tielens & Hagen 1980). Upon heating from the 
accreting central protostar, the grain ices sublimate inject¬ 
ing those molecules in the gas-phase and starting in this 
way what was then called a warm gas-phase chemistry. This 
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paradigm has lasted for almost two decades, until new ob¬ 
servations (e.g. Cazaux et al. 2003; Bottinelli et al. 2007; 
Requena-Torres et al. 2006), laboratory experiments (Gep- 
pert et al. 2007) and theoretical studies (Horn et al 2004) 
challenged it. The warm gas phase chemistry paradigm left 
the place to the grain surface chemistry. 

In the new paradigm, grain surface chemistry is not 
only responsible for the hydrogenated molecules of the pre¬ 
collapse phase, but also for the (almost) whole set of ob¬ 
served COMs. The basic idea of this class of models (Garrod 
& Herbst 2006) is that radicals trapped in the iced mantles 
acquire mobility and react forming COMs when the dust 
temperature reaches ~ 30 K. The origin and quantity of 
the trapped radicals depends on the adopted model. They 
might be the pieces of iced-species (e.g. methanol) broken 
by the FUV photons produced by the interaction of cos¬ 
mic rays with H 2 (e.g. Garrod & Herbst 2006; Chang & 
Herbst 2014), or, alternatively, they are the result of the in¬ 
complete hydrogenation of the simple mother species (e.g. 
Taquet et al. 2012). Recent observations, however, show the 
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Figure 1. Scheme of the added reactions that form DME, MF 
and methoxy (boxes with grey large contours). 

presence of some COMs -notably methoxy (CH 3 O), methyl 
formate (HCOOCH 3 , MF) and dimethyl ether (CH 3 OCH 3 , 
DME)- in regions where the dust temperature is less than 
30 K: pre-stellar cores (Bacmann et al. 2012; Vastel et al. 
2014) and cold envelopes of low mass protostars (Oberg et 
al. 2010; Cernicharo et al. 2012; Jaber et al. 2014). Two re¬ 
cent models have been published trying to explain these new 
observations, assuming that the detected COMs are present 
inside the cold and dense core (Vasyunin & Herbst 2013, 
hereinafter VH2013; Chang & Herbst 2014). 

More recent observations towards the pre-stellar core 
L1544 (Bizzocchi et al. 2014; Vastel et al. 2014) added a 
crucial information: the methanol emission originates from 
the outer shell of the cold core, where the H 2 density is 
~ 3 X 10'* cm“® and the temperature is ~ 10 K, and not 
from the inner and denser region, as previously assumed. 
Therefore, methanol is ejected from the grain mantles into 
the gas because of some non-thermal ice desorption, as it 
was previously claimed for water (Caselli et al. 2012). Fi¬ 
nally, Vastel et al. (2014) argued that the COMs emission 
detected in LI544 also comes from the same outer shell and 
speculated that this is likely the case also for the other cold 
sources where COMs have been detected (Bacmann et al. 
2012; Oberg et al. 2010; Cernicharo et al. 2012). 

Following these last findings and building up on the 
VH2013 model, we propose here a new model where gas 
phase reactions, triggered by the non-thermal desorption of 
methanol from the ices, play a major role in the formation 
of methoxy, MF and DME in cold regions. Specifically, we 
consider a set of gas phase reactions, which were previously 
overlooked or whose coefficients we have refined based on 
recent experimental and theoretical works. 


2 MODEL AND ADOPTED CHEMICAL 
NETWORK 

2.1 The model 

We used the GRAINOBLE code (Taquet et al. 2012, 2014) 
to compute the mantle species from hydrogenation of CO 


and O. Briefly, GRAINOBLE couples the gas-phase and 
grain-surface chemistry with the rate equations approach 
introduced by Hasegawa et al. (1992), and takes into ac¬ 
count the multilayer and porous nature of the ice mantles. 
Species accrete on the grain surfaces and are desorbed by 
cosmic ray-induced heating of grains (Hasegawa & Herbst 
1993) and by a non-thermal process that injects into the 
gas phase 1 % of the species formed on the surface, follow¬ 
ing a process called “chemical desorption” (Garrod et al. 
2007). We followed the chemical composition evolution of a 
gas with a H density of 6 x 10* cm“® and a temperature 
of 10 K (to simulate the conditions in L1544; Vastel et al. 
2014). We started with a gas having the initial abundances 
quoted by Wakelam & Herbst (2008) for the diffuse clouds 
(their Table 1, column EA2; note that hydrogen is initially 
molecular). We follow the chemical evolution for 10® yr. 

In addition to this “standard” model, we run a few more 
models to evaluate the impact on the results when some cru¬ 
cial parameters are changed. Specifically, we considered the 
following cases: (i) a 0 % or 10 % non-thermal desorption frac¬ 
tion, the latter chosen based on the experiment by Dulieu et 
al. (2013), although new experiments by Minissale & Dulieu 
(2014) suggest a dependence of this value on the substrate 
and other parameters; (ii) a gas with a 5 times larger den¬ 
sity, and with 20 K temperature; (iii) a high and low Cl 
and F abundances to evaluate their impact on the results, 
as some important reactions introduced in this work depend 
on them (see next §); (iv) the rate coefficient of reaction ( 2 ) 
of Table 1 reduced by a factor 10 (see below). 

2.2 The chemical network 

The surface chemical network is based on the one presented 
by Taquet et al. (2013) but does not include the formation 
of deuterated species. We used the binding energies listed 
in Taquet et al. (2014) and the energy barrier for the CO 
and H 2 CO hydrogenation recently computed by Rimola et 
al. (2014). 

For the gas phase chemistry, we used the 
OSU2009 network (Harada & Herbst 2008; 
http://faculty.virginia.edu/ericherb) as a basis, and 
added a set of new reactions for the formation of MF, DME 
and methoxy, and involving atomic chlorine and fluorine 
for the conversion of DME into its radical CH 3 OCH 2 . The 
scheme of the chemical reactions added in this work is 
shown in Fig. 1 and their list is reported in Table 1. We 
comment the novelties of the adopted network, for each 
species, in the following. 

Methoxy is mostly formed by reaction 1. This reac¬ 
tion was also considered by VH2013, with a rate coef¬ 
ficient equal to 4 x 10“** cm®s“*. We instead adopted 
the value calculated by Shannon et al. (2013) at 20 
K, 3 X 10 *** cm®s *, which, according to their sugges¬ 
tion, is likely a lower limit. Note that the reactions in¬ 
volving the methoxy reported by VH2013 in Table 1 
are not in line with the cited NIST Chemical Kinet¬ 
ics Database {http://kinetics.nist.gov/kinetics/Search.jsp), 
where the rate coefficients are different from those quoted 
by VH 20 lQ 

^ The rate coefficients reported by VH2013 for the reactions 
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Reaction 


Rate coefficient References Sz 
(cm^ at 10 K) Notes 


(1) OH + CH 3 OH-5- CH 3 O + H 2 O 3.0x10-“ 

(2) CH 3 O + CH 3 -s. CH 3 OCH 3 + photon 3.0x10-“ 

(3) O + CH 3 OCH 2 -s> HCOOCH 3 + H 2.0x10-“ 

(4) F + CH 3 OCH 3 -5- CH 3 OCH 2 + HF 2.0x10-“ 

(5) Cl + CH 3 OCH 3 -)■ CH 3 OCH 2 + HCl 2.0x10-“ 


1 

see text 
2, 3 
2 
4 


Table 1. List of the added reactions that form DME, MF and methoxy. References: 1- Shannon et al. 2013; 2- Hoyermann & Nacke 
1996; 3- Song et al. 2005; 4- Wallington et al. 1988. 


DME is mainly formed by reaction 2. The rate coefficient 
for the radiative association of methyl and methoxy is not 
kwown. Recent theoretical work by Vuitton et al. (2012) 
suggests that radiative association can be very efficient 
when two relatively large radicals interact and form a stable 
molecule, if two-body exothermic channels are not available 
to the system. VH2013 quote the presence of a competitive 
exothermic channel leading to H 2 CO -I- CH 4 and, therefore, 
employed a rate coefficient at 10 K of 2.6 x lO-^'^ cm^s-^. 
However, a recent study by Sivaramakrishnan et al. (2011) 
has clearly pointed out that the reaction CH 3 O -|- CH 3 —>■ 
H 2 CO + CH 4 is a direct abstraction process, not correlating 
with the bound intermediate CH 3 OCH 3 . For this reason, we 
have increased the rate coefficient for reaction 2 by an order 
of magnitude, as there are no alternative two-body reaction 
channels and the hot DME formed after the association of 
CH 3 and CH 3 O can only dissociate back or stabilize by ra¬ 
diative association. Yet, we have run also a model with the 
reaction 2 rate coefficient a factor ten smaller to evaluate 
the impact on the results. 

MF is formed by reaction 3, namely the oxidation of 
CH 3 OCH 2 . This radical is formed by several reactions in¬ 
volving DME (Fig. 1) and common atomic and radical 
species, such as atomic chlorine and fluorine following re¬ 
actions 4 and 5. The employed rate coefficients for reactions 
3, 4 and 5 have been determined by Hoyermann & Nacke 
(1996) and Wallington et al. (1988) at 298 K. Reactions 3 
to 5, however, do not have an entrance barrier and we can 
retain their value also at much lower temperatures (as a mat¬ 
ter of fact, a moderate increase of the rate coefficient with 
decreasing temperature is expected in this case). These re¬ 
actions, well characterised in laboratory experiments, were 
not included in previous astrochemical models. We empha¬ 
sise that they provide a direct link between DME and MF, 
being the former a precursor of the latter. 

In addition, the standard loss mechanisms of COMs 
present in OSU (such as the reactions with H*^, He^ etc) 
have all been retained and we have added similar loss path¬ 
ways (with similar rate coefficients) also for the new species, 
such as CH3OCH2, which are not present in OSU. In ad¬ 
dition, we considered the standard reactions of atomic F 
(Neufeld et al. 2005) and employed the recent determina¬ 
tion of the rate coefficients for F -f H 2 by Tizniti et al. 
(2014). As for the reactions of atomic Cl, we note that the 

CH 3 O + H and CH 3 O -I- O are erroneous: the NIST values 
are 3 X IQ-^^ and 2.5 X IQ-cm^s-^, respectively, rather than 
1.6 X IQ-^® and 1.0 X IQ-In addition, the dominant reaction 
channel for CH 3 O + H is the one leading to H 2 CO + H 2 and not 
the one leading to CH 3 -|- OH, as claimed by VH2013. 


Cl -I- H 2 reaction is endothermic (Balucani et al. 2004, and 
references therein) and cannot occur in the conditions of the 
ISM (Neufeld et al. 2009). Other reactions between chlorine 
and small saturated molecules are known to be very inef¬ 
ficient. Therefore, Cl is left to react with larger molecules, 
such as those of interest here. Please note that the Cl role in 
the COMs chemistry has been overlooked so far, although 
it can (only) react with large molecules. 


3 RESULTS 

Figure shows the evolution with time of the abundances 
of methoxy, DME and MF, plus CH3 and methanol. In the 
standard case, the methoxy, DME and MF gaseous abun¬ 
dances reach a peak of ~ 2 x ~ 2 x and 

~ 5 X respectively, at ~ 10® yr and then decrease 

to lower values because they freeze-out onto the grain man¬ 
tles. The peak occurs at ~ 10® yr because of the formation 
of methanol on the grains, which occurs when CO starts to 
freeze-out onto the grain mantles. A fraction of the methanol 
is then released in the gas phase and this triggers the re¬ 
actions of Fig. 1, which form methoxy, DME and MF. In 
the same panel, we also show the predicted abundances ob¬ 
tained without considering the new reactions of Table 1. As 
expected, the predicted abundances of DME and MF are 
respectively about three and one order of magnitude lower. 

The panel b of Fig. shows the impact on the abun¬ 
dances when the chemical desorption is varied. The increase 
by a factor 10 causes a proportional increase of gaseous 
methanol and methoxy, but a larger increase of DME and 
MF abundances, because methoxy is the bottleneck to the 
formation of DME. Caused by the dominance of reactions 

4 and 5, the MF abundance increases by the same factor 
than DME. Conversely, putting the chemical desorption to 
zero kills the Fig. 1 gas phase reactions as standard cosmic- 
ray induced desorption is not effective to desorb methanol 
from the ices (as GRAINOBLE considers layer and not bulk 
chemistry in mantles). 

The increase of density by a factor of 5, shown in in 
the panel c of Fig. causes a shorter time of residence 
of the molecules in the gas phase but only slightly affects 
the absolute abundance of the species. On the contrary, the 
increase of the temperature to 20 K causes a decrease of the 
methoxy, DME and MF peak abundances. This is due to 
the higher abundance of gaseous atomic H (whose density 
is about 1 cm-® at 10 K). At 20 K, the density of atomic 
H increases by one order of magnitude due to the higher 
evaporation rate. 

The decrease by a factor 10 of the rate coefficient of 
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Figure 2. Predicted abundances (with respect to H nuclei) of MF (HCOOCH 3 ), DME (CH 3 OCH 3 ) and methoxy (CH 3 O) as a function 
of the time, plus some species that intervene in their formation, CH 3 and methanol (CH 3 OH). The upper left panel (a) refers to the 
standard case (dashed lines) and the case without the new reactions of Table 1 (dotted lines). The upper right panel (b) is the same than 
(a) but with a chemical desorption rate equal to 0 (dotted lines) and 10% (dashed lines) respectively. The lower left panel (c) shows a 
gas with a H density of 3 X 10® cm“® (dotted lines) and temperature 20 K (dashed lines). The lower right panel (d) shows the standard 
case but with the rate of reaction 2 of Table 1 10 times lower (dotted lines) and with low (see text) Cl and F abundances (dashed lines). 


reaction 2 leads to a similar decrease of the DME abundance. 
It also leads to the MF abundance forming earlier (because 
of other reactions than reaction 3 of Table 1) than DME and 
reaching a similar low abundance. Finally, decreasing the Cl 
and F elemental abundance (to the low values quoted by 
Graedel et al. 1982) results in a similar behavior of the MF 
abundance, as reactions 4 and 5 of Table 1 do not form 
enough CH 3 OCH 2 . Note that F is quickly (~ 10® yr) locked 
into HF, while Cl remains mostly atomic until the freeze-out 
kicks on (~ 10® yr). Therefore, reaction 5 of Table 1 plays 
a major role than reaction 4 in the formation of CH 3 OCH 2 
and, consequently, in the formation of MF from DME. 


4 DISCUSSION 

Whether the proposed new model reproduces the COMs 
abundances observed in cold objects depends on the inter¬ 
pretation of where their emission comes from. Vastel et al. 
(2014) and Bizzocchi et al. (2014) showed that methanol 
emission originates in a shell at the border of the condensa¬ 
tion of the pre-stellar core L1544. These authors claim that 
methanol is non-thermally desorbed from the grain mantles 


and that, given the lower gas density, the rate of desorption 
overrun the rate of freezing-out. Since the same argument in 
principle applies also to the COMs, Vastel et al. (2014) sug¬ 
gested that COM emission originates from the same outer 
shell. If this is true, the derived abundance of the COMs 
detected by Vastel et al. (2014) would be 20 times larger 
than if they originated in the cold and dense inner part of 
the L1544. Evidently, a similar factor may affect the COMs 
abundance derived in the other cold regions (Bacmann et 
al. 2012; Cernicharo et al. 2012). 

Vastel et al. (2014) estimated a gaseous methanol abun¬ 
dance in the FUV-illuminated shell of ~ 6 x 10“®, and an 
upper limit to the methoxy, DME and MF abundances of 
~ 2 X 10~®®, ~ 2 X 10~®® and ~ 2 x 10“®, respectively. 
These numbers agree very well with the predictions of our 
standard model. On the contrary, the 10% chemical desorp¬ 
tion model predicts too large abundances. Unfortunately, 
the upper limit on the MF does not allow to put constraint 
on the role of Cl and F, nor on the value of the reaction 2 rate 
coefficient. Although we cannot directly compare our predic¬ 
tions with the absolute abundances quoted by Bacmann et 
al. (2012) and Cernicharo et al. (2012) because the assumed 
H 2 column densities may be overestimated, the DME over 
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MF abundance ratios are compatible with our predictions. 
On the contrary, the methoxy abundance measured by Cer- 
nicharo et al. (2012) is 5 times lower than the DME and MF 
abundances, at odds with our predictions. 

Constraints on the reactions 3 to 5 come from mea¬ 
sures of the DME and MF abundances in different objects. 
Specifically, Jaber et al. (2014) showed that there is a linear 
correlation between these two quantities over a six orders 
of magnitude range, confirming a direct link between DME 
and MF. The reactions 3 to 5 that we propose here explain 
in a natural way that link. Similarly, we predict the methoxy 
to be correlated with the methanol abundance, as it is the 
product of the reaction of CH 3 OH with OH. This reaction 
has been recently studied and found to be much more ef¬ 
ficient than previously expected (Shannon et al. 2013). Fi¬ 
nally, we predict a possible correlation between methoxy and 
DME, although the correlation could be broken if CH 3 be¬ 
comes the bottleneck of the reaction 2. This correlation is 
based on the associative reaction between CH 3 and CH 3 O, 
a likely process at 10 K, which, however, has never been in¬ 
vestigated (if- Since the determination of the product yield 
and rate coefficients for this reaction will be extremely dif¬ 
ficult, if not impossible, in laboratory experiments, we urge 
our colleagues to carry out theoretical calculations of the 
same kind as those carried out by Vuitton et al. (2012) for 
other species. 


5 CONCLUSIONS 

We propose a new model to explain the formation of MF, 
DME and methoxy, based on recent experimental and the¬ 
oretical works. The major conclusions of this work are: 

i) the non-thermal desorption of iced methanol starts a series 
of gas phase reactions leading to the formation of methoxy, 
DME and ME; 

ii) Methoxy is formed via the methanol-hydroxyl reaction; 

iii) DME is mainly formed via the radiative association re¬ 
action of methoxy with CH 3 ; 

iv) MF is formed from DME via the oxidation of CH 3 OCH 2 , 
a reaction overlooked by previously models; 

v) MF is, therefore, the daughter of DME, which explains 
the observed correlation (e.g. Jaber et al. 2014). 

In conclusion, grain surface chemistry certainly plays a 
role, for example in forming hydrogenated species during the 
prestellar phase (Tielens & Hagens 1980; Caselli & Ceccarelli 
2012), but not necessarily in the formation of all COMs. 
Gas phase and grain surface chemistry will have to share 
the reign of the formation of COMs. 
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